Random phases in Bose-Einstein condensates with higher order nonlinearities 
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We present a statistical description of Bose-Einstein condensates with general higher order non- 
linearities. In particular, we investigate the case of cubic-quintic nonlinearities, of particular interest 
for dilute condensates. The implication of decoherence for the stability properties of the condensate 
is discussed. 
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The basic concept of macroscopic quantum states, 
such as Bose-Einstein condensates (BECs)^ and lately 
also Fermi condensates^, has caught the interest of the 
physics community, both due to the nature of the con- 
cept itself but also since the possibilities to perform new 
and exciting experiments was early recognized 1 . BECs 
are normally described by the Gross-Pitaevskii (GP) 
equation^, in which the cubic nonlinearity represents two- 
body forces between the bosons in the condensate. There 
are numerous works on the theoretical foundations and 
implications o&££. When the scattering length is posi- 
tive, the possibility of dark solitons is given via the GP 
equation. Such dark solitons are stable towards pertur- 
bations in one dimension, while the multi-dimensional 
case is more complex. Such dark solitons have also been 
found experimentally 7 . In the case of a negative effective 
scattering length, the GP equation admits bright soli- 
tons, which are prone to collapse in dimensions larger 
than one. Such bright solitons have also been experi- 
mentally obtained, both in terms of trains of solitons 8 ., 
as well as single solitary structures^. The dynamics of 
these bright soliton structures has also been analyzed 
both analytically 10 and numerically 11 . 

If three-body interactions are taken into account, 
higher order nonlinearities will modify the GP equation. 
We may write this as a generalized NLSE of the formic 



ihdtip+ — V 2 ip + a\ip\ 2 ip + /3M 4 ip = 0, (1) 
2m 

where a and (3 in general is complex-valued, and ip is 
the condensate wave function. As goes to zero, we 
regain the GP equation. The real parts of a and (3 corre- 
sponds to elastic collisions within the condensates, while 
the imaginary part appears due to inelastic scattering 13 . 
In what follows we will neglect the collisional losses, and 
assume that a and (3 is real, and can take on both pos- 
itive and negative values. In fact, the coefficient a is 
proportional to a, the scattering length, which can be 
tuned to take on both negative and positive values (see, 
e.g. Refs. Q and ll7^ . while (3 is proportional to a 4 for 
dilute systems 15 : 16 ! 17 . Thus, in BECs where a is large, 
the quintic contribution to Eq. Q may become signifi- 
cant. Equation Q can be demonstrated to have solitary 



solutions in one and two dimensions^, and appears not 
only in the physics of BECs, but also in, e.g. nonlinear 
optics^. Furthermore, Eq. QJ is a special case of the 
equation 

h 2 

iWt^ + —W 2 ip + U(\'>p\ 2 )ip = 0, (2) 
2m 

where U in general is a complex-valued function of the 
norm of the wave function squared. 

The stability of solutions to Eq. (JJJ towards coher- 
ent perturbations was analyzed by Shukla and Yu**, 
where the growth rate for the modulational instability 
was found. This growth rate signifies the onset of purely 
growing perturbations, and is thus an important indi- 
cation of the possibility of solitary solutions. However, 
the effects of incoherence, e.g. a random phase in the 
wave function, may significantly alter the modulational 
instability and therefore also the onset of inhomogeneity 
growth, and is an important issue (for a discussion, see, 
e.g. Refs. H3 and IH. A very direct approach in analyz- 
ing the effects of partial coherence lies in the Wigner 
formalism 2 ^ 2 ^. This approach has found uses in the 
study of surface gravity waves in fluids* 2 * and electro- 
magnetic waves in nonlinear medial, in quantum statis- 
tical mechanics^ 6 -, in nonlinear optics 2 ^, and in quantum 
plasmas 28 . Here we will apply the Wigner formalism to 
the problem of partial coherence in the modulational in- 
stability of higher order nonlinear BECs. 

In order to analyze the statistical properties of Eq. QJ , 
we may introduce the Fourier transform of the two-point 
correlation function of the wave function, i.e. the Wigner 
function, according to 



F(t,r,p) 



(27T/l) f 



ip-i/h 



(^*(r+£/2,f)^(r-£/2,i)) 
(3) 



for the wave function ip, where s denotes the dimension- 
ality of the problem at hand, the asterisk is the complex 
conjugate operation, and the angular bracket denotes the 
ensemble average. The Wigner function corresponds to a 
generalized distribution function for the bosons, and by 
applying the time derivative to Eq. JSJ and using Eq. J2J), 



2 



one finds the Vlasov-like equation 

d t F + -p.VF+ lu(\tp\ 2 ) sin f J V ■ V P ] F = 0. (4) 
to h \Z J 

where the sin-operator is defined in terms of its Taylor 
expansion and the arrows denote the direction of opera- 
tion. In the case of of a cubic-quintic nonlinearity, such 
as in Eq. (JJJ, we obtain 



t F + -p-VF+l(a\i;\ 2 +P\ij\ i )sm fj V ■ V P )f = 
to n \z J 

(5) 

Moreover, the modulus square of the wave function is 
given by 



= J dpF(t,r,p). 



(6) 



For the sake of clarity, we now focus on the one- 
dimensional case. The stability of equation J^J) can 
be analyzed using a linearization procedure. Letting 
F = Fq(p) + f(p)exp(ikz — icvt), where |/| <C F , we 
linearize Eq. in order to obtain the nonlinear disper- 
sion relation 



1 = - 



to dU 
hk d\ipo\' 



dp 



F {p + hk/2)-F {p-hk/2) 
p — tom/k 



Equation Q is the general dispersion relation for matter 
waves taking into account higher order nonlinearities. 

For the case of a cubic-quintic nonlinearity, Eq. Q 
reduces to 

F (p + hk/2)-F a {p-hk/2) 



1 = - — (a+2/3h/> | 2 ) / dp- 



p — Lum/k 



For the monochromatic wave case, i.e. Fq(p) 
IV'ol^tj' — Po), we obtain the dispersion relation^ 



CO 



p Q k 



± 



f h 2 k 4 
4m 2 



k 2 M 



■(a + 2f3\ip \ 



1/2 



(9) 



from Eq. JHJ). Setting (3 = in the dispersion relation 
ijffil. we obtain the Bogolubov expression^ for the ele- 
mentary excitations of the BEC. The standard method 
employed in obtaining the result 10 with (3 = is to 
set ip = ipo + ipi , where \ipi\ <^ |^>o| and "00 is the back- 
ground state, and linearizing Eq. (Tjl. after which the 
equation may be split into its real and imaginary part 
and harmonically decomposed. We note that the Wigner 
approach presented in this paper is equivalent in to the 
Bogolubov method in the monochromatic limit. Letting 
lo = pok/m + in the above equation, we obtain the 
modulational instability growth ratei^ 



•(a + 2/3|Vo| 2 )- 



h 2 k 4 

Am 2 



1/2 



(10) 



If the waves are not exactly monochromatic, but have 
a spectral broadening due to, e.g. a random phase in the 
background wave function, we may model equilibrium 
condensate spectrum by a Lorentzian distribution 30 



F (p) 



PT 



(P 



Po) 2 +P 2 T 



(11) 



where pt denotes the width of the distribution. The 
Lorentzian distribution solves Eq. JSJ , as well as the more 
general (0J, and is thus a valid perturbation background. 
We note that the phase fluctuations may stem from a 
variety of perturbations, e.g. thermal effects or quantum 
fluctuations^. The dispersion relation J7J) for this case 
then turns out to be 



p k 



± 



h 2 k 4 



dU 



-,1/2 



4to 2 to dlV'ol 2 
which gives the growth rate 
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k 2 dU 
to d|i/>o| ; 



h 2 k 4 

4^ 



i^, (12) 

TO 



1/2 



Prk 

TO 



(13) 



The growth rate l|13|) is valid for a general nonlinear- 
ity. We note that if dU/d\ip \ 2 < 0, there is no mod- 
ulational instability growth, and the perturbations are 
damped. Thus, a minimum requirement for a positive 
growth rate for an arbitrary nonlinearity in Eq. @ is 
that dU/d\tp \ 2 > 0. 

In the case of a cubic-quintic nonlinearity C/(|?/;| 2 ) = 
\ip\ 2 {a + P\ip\ 2 ), Eq. U2} reduces to 



Pok 

Ul = ± 

TO 



(8) uj = — 



h 2 k 4 fcW , |2 .1 1/2 . PT k 
-7—2 + -1 , 

4m z to to 

(14) 

which gives a purely growing modulational instability 
whose growth rate is 



7 : 



k 2 \U 



m 



ia + 2p\?p \ 2 ) 



h 2 k 4 



4m 2 



1/2 



^ (15) 

TO 



Comparing with Eq. I|15l) with Eq. IjlUI) , one clearly sees 
the damping character of the spectral broadening term. 
We note that since a and (3 may be positive or negative, 
independent of each other, the instability properties cru- 
cially depends on the nonlinear terms in the expression 
(|13|l . If a, (3 < (the defocusing case), modulational 
instability growth is not possible, and this corresponds 
to the well-known stability of dark solitary solutions. 
However, if we have a < (defocusing cubic nonlinear- 
ity), while (3 > (focusing quintic nonlinearity), we may 
have a new instability regions, not present in the Gross- 
Pitaevskii equation. If a > (focusing cubic nonlinear- 
ity) and (3 < (defocusing quintic nonlinearity) we will 



3 



have a damping in the growth rate due to three-body in- 
teraction. Note that this damping is quite different from 
the dissipation due to inelastic three-body scattering, and 
has more the character of Landau dampingSi. It is also 
clear that the case a,(3 > (the focusing case) gives the 
maximum instability growth rate. However, even in this 
case, a very broad spectral distribution of the BEC, due 
to e.g. thermal noise, may quench the growth rate con- 
siderably, even removing it all together. We would like 
to stress that the Wigner method presented in this paper 
is a very general approach to partial coherence and spec- 



tral broadening. Thus, the method is appropriate for a 
variety of spectral distributions, e.g. Gaussians, as well 
as the Lorentizian , but this will be pursued in future 
research. 

To summarize, we have presented a perturbation anal- 
ysis of the statistical properties of the generalized Gross- 
Pitaevskii equation by using the Wigner formalism. 
In the case of a random phase of the background wave 
function of the condensate, we find that the spectral 
broadening gives rise to a reduced growth rate, as com- 
pared to the mono-chromatic case. 
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